Mott and Wigner-Mott transitions in doped correlated electron systems: effects of 
superlattice potential and inter-site correlation 
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We introduce the notion of superstructure Mottness to describe the Mott and Wigner-Mott tran- 
sition in doped strongly correlated electron systems at commensurate filling fractions away from 
one electron per site. We show that superstructure Mottness emerges in an inhomogeneous electron 
system when the superstructure contains an odd number of electrons per supercell. We argue that 
superstructure Mottness exists even in the absence of translation symmetry breaking by a super- 
lattice, provided that the extended or intersite Coulomb interaction is strong. In the latter case, 
superstructure Mottness offers a unifying framework for the Mott and Wigner physics and a non- 
perturbative, strong coupling description of the Wigner-Mott transition. We support our proposal 
by studying a minimal single-band ionic Hubbard t-U-V-A model with nearest neighbor Coulomb 
repulsion V and a two-sublattice ionic potential A. The model is mapped onto a Hubbard model 
with two effective "orbitals" representing the two sites within the supercell, the intra and interorbital 
Coulomb repulsion U and U' ~ V, and a crystal field splitting A. Charge order on the original lattice 
corresponds to orbital order. Developing a cluster Gutzwiller approximation, we study the effects 
and the interplay between V and A on the Mott and Wigner-Mott transitions at quarter-filling. 
We provide the mechanism by which the superlattice potential enhances the correlation effects and 
the tendency towards local moment formation, construct and elucidate the phase diagram in the 
unifying framework of superstructure Mottness. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.30.+h 



I. INTRODUCTION 

The Mott physics, i.e. the transition from band to 
localized states due to interaction, plays an important 
role in materials containing strongly correlated electrons. 
The usual condition for Mott transition on a lattice re- 
quires half-filling, i.e. one electron per site. Due to 
Coulomb blocking, the electrons are localized and be- 
have as local moments when the on-site Coulomb repul- 
sion U becomes larger than the bandwidth of the kinetic 
energy. The transition metal oxides and rare earths and 
actinides, involving narrow bands from the atomic d and 
/ orbitals, arc known to exhibit Mottness - a remarkable 
set of electronic properties due to the proximity to Mott 
transitionii^. 

Thus far the study of Mott transition and the associ- 
ated Mottness has mostly focused on uniform systems. In 
this paper, we propose and develop the notion of super- 
structure Mottness to describe the Mott physics for inho- 
mogeneous electronic states. This is important because 
strongly correlated electron systems have a propensity to- 
wards inhomogeneous phases as a result of the frustrated 
kinetic energy. Examples include the transition metal 
chalcogenidcs that form charge density waves (CDW)^i^ 
and the high-Tc cuprates that exhibit stripe^ and checker- 
board orders^. The superstructure in strongly corre- 
lated systems can also come from structural distortions 
as in Ca-substituted strontium ruthenates-^ or chemical 
dopants as in the sodium cobaltates^. 

Let us consider a commensurate superstructure in two 
dimensions with a, p x q supercell. An electron density 
n* — £/pq corresponds to i electrons per supercell. If £ is 
an odd integer, the condition for the superstructure Mott 



transition is satisfied by an odd number of electrons per 
supercell. More precisely, ior £ = 2k + l, k of the M = pq 
subbands are completely filled, while the (fc + l)-th sub- 
band, hosting one electron, is at half-filling and can un- 
dergo Mott transition. Note that partial filling of the 
subbands leading to the semi-metallic states, which usu- 
ally happens close to integral filling near the top and bot- 
tom of the subbands for an even number of electrons per 
supercell, will not occur here for a generic superlattice 
potential. As a result, superstructure Mottness allows 
the Mott transition to occur far away from one electron 
per site. For example, electron correlation is expected to 
be weak in the sodium rich phases of Naa.Co02 because 
of the proximity to a band insulator at x = 1. To the 
contrary, strong correlation effects are observed, possibly 
due to the superstructures induced by Na dopant o^'^'^°i^^ 
that allow superstructure Mottness. Another example is 
the recent observation of orbital selective Mott transition 
in Ca2-2:Sr2:Ru04 near x = 0.3^^; the dxy orbital hosting 
1.5 electrons becomes Mott localized due to the v^ x ^/2 
superstructure induced by RuOg octahedral distortion^ 
such that there are 3 electrons per supercell. A third 
example is that the vTS x vTS CDW reconstruction on 
the surface of lT-TaSe2 leads to one electron per super- 
cell and causes a Mott transition at low temperatures 
which has been observed recently by scanning tunneling 
spectroscopjii. 

In addition to Mott localization, another paradigm 
in the physics of strong electronic correlation is Wigner 
crystallization which, in contrast, is due to strong long- 
range Coulomb interaction V that causes electrons to 
crystallize (charge order) by Coulomb jamming. The 
electron crystallization into a charge ordered Wigner lat- 



(a) A 




(b) A ^ 



A band 



£band 




FIG. 1: (Color online) (a) and (b) Schematics of band evolu- 
tion with interaction. 



tice is another form of electronic superstructures. Wigner 
localization happens even classically for low density elec- 
trons in the continuum. What happens to lattice elec- 
trons when both U and V are large and away from half- 
filling, i.e. the Wigner-Mott transition, is largely unex- 
plored theoretically because of the lack of nonperturba- 
tive many-body methods for treating the finite-range, 
inter-site Coulomb repulsion V . The basic difficulty is 
that the Hilbert space becomes nonlocal in the pres- 
ence of the inter-site V, which defies the strong-coupling 
approaches based on the single-site Hilbert space con- 
taining four states representing an empty site, an up- 
spin electron, a down-spin electron, and a doubly occu- 
pied site. Previous approaches used a single-site Hilbert 
space Gutzwiller projection or dynamical mean field 
theory (DMFT) to treat the onsite C/, but decoupled 
the inter-site y-term using the weak-coupling Hartree 
approximatio n^^'^ 1^^ which becomes unphysical when V 
is large since the electrons would avoid the high self- 
energy cost by not sitting close together. 

It turns out that the notion of superstructure Mottness 
encompasses the physics of Wigner-Mott transition in a 
unified framework. Our basic insight is to envision the 
lattice as being covered by clusters or supercells and think 
of the latter as artificial atoms with a number of effective 
orbitals representing the sites contained in each cluster. 
This maps the problem to that of a multi-orbital Mott- 
Hubbard system where U and V play the role of intra- 
and inter-orbital interactions. Charge ordered states on 
the original lattice simply correspond to orbital order in 
the artificial atoms. This remarkable mapping, although 
approximate, allows us to treat the effects of V using the 
nonperturbative strong coupling approach for the multi- 
orbital Hubbard model and thus one is led to embrace 
a possible emergent Mott transition with increasing U 
and V at the filling fraction of one electron per cluster 
(supcrcell). Such a phenomenon of localization and lo- 
cal moment formation into clustered superstructures is 
another form of superstructure Mottness that unifies the 
Mott and Wigner physics. 

The above picture amounts to truncating the nonlocal 
Hilbert space as a product of that of the clusters or super- 



cells. The interactions between the electrons in the neigh- 
boring clusters are treated on equal footing. We extend 
the well established Gutzwiller projected wavefunction 
approach and develop a cluster Gutzwiller method. Ex- 
plicit model calculations are performed to test the ideas 
of superstructure Mottness. Specifically, we study a min- 
imal one-band, ionic Hubbard model with intersite cor- 
relations, i.e. the t-U-V-lS. model. In addition to the on- 
site Hubbard U and the nearest neighbor (NN) Coulomb 
interaction V , an onsite two-sublattice potential A^ is in- 
cluded for possible external potentials due to lattice dis- 
tortion and dopant ions. We map the model to a generic 
two-orbital Hubbard model on a superlattice that con- 
tains two sites per supercell. The Hubbard U serves as 
the intra-orbital Coulomb U\ the inter-site V becomes 
the inter-orbital repulsion U' ~ V, whereas the superlat- 
tice potential A^ turns into an effective crystal field Aq, 
with a — A,B labeling the two "atomic" orbitals with 
a crystal field splitting A = A^ — A^. Using the clus- 
ter Gutzwiller projection of the multioccupancy states 
in a supercell, in analogy to the multiorbital Gutzwiller 
projection methodic, we study the ground state phase 
diagram at quarter filling in the parameter space of U, 
V, and A. Fig. [1] shows the schematics of two classes 
of localization transitions with increasing U and V . For 
A = 0, the two bands remain degenerate and undergo 
the transition simultaneously from a uniform correlated 
metal to a Mott insulator without charge (i.e. orbital) 
ordering. For A 7^ 0, we find a transition from a CDW 
metal involving both bands to a charge ordered insulator, 
where one of the two bands undergoes an orbital selec- 
tive Mott transitionii. We elucidate the mechanism by 
which the superlattice potential enhances the correlation 
effects and the tendency towards local moment forma- 
tion, and reveal a deeper connection among the strongly 
correlated inhomogeneous electronic states, the Wigner- 
Mott physics, and the multiorbital Mott physics that can 
all be united under the notion of superstructure Mot- 
tness. 



The rest of the paper is organized as follows. In sec- 
tion II, we discuss the model and the mapping to an 
equivalent multiorbital Hubbard model on the super- 
structure. In section III, we describe the virtual cluster 
Gutzwiller approximation (VCGA) and study the super- 
structure Mott transition in the ionic Hubbard model, 
i.e. the t-U-A model at quarter-filling, in the absence of 
intcr-site Coulomb correlation V, In section IV, we ex- 
tend the VCGA to the case of finite intcr-site correlation 
V and study the Wigner-Mott transition in the absence 
of the ionic potential A, i.e. the t-U-V model at quater- 
filling, using the ideas of superstructure Mottness. The 
results will be compared to the weak-coupling Hartree 
approximation. The general case of finite A and V is 
presented in section V followed by a brief summary in 
section VI. 



II. MODEL AND HEURISTIC DISCUSSION 



The Hamiltonian of the t-U-V-A model is given by 



H 



i,3 



{Uj - Ai6^j)c\^Cj„ + uY^ ^iT^a - 



'V'^fiifij, 



(1) 

where c]^ creates an electron at site i of spin ct on a 
square lattice. We consider electron hoppings between 

/uj)) ~ t' . Repeated spin 



NN t 



{*j> 



t and next NN t 



indices are summed. Eq. ([T]) is an extended version of 
the ionic Hubbard model with NN intersite repulsion V 
and a lattice ionic potential A^. We focus on the simplest 
case where A^, when present, has a superstructure with 
two-sites per supercell. 

Consider the lattice as a complete coverage by a two- 
sublattice {A and B) superstructure. The two sites per 
supercell (cluster) can be viewed as two effective or- 
bitals. The hopping terms can be written in the equiv- 
alent form of two bands having dispersions ^aa(fc) = 
— 2f'(cosfca;-|-cosA:y) with an inter-band (a ^ /3) hopping 
^Q,/3 = — 4icos(A:a;/2)cos(fcy/2), where a,(} = A,B label 
the two orbitals and k £ the reduced Brillouin zone of 
the reciprocal lattice. The Hamiltonian can be rewritten 
in a suggestive form 



Q,/3,fe aj 

+ ^ Aaflj^a + U''^ ni,Ani,B 
aJ I 

+ V^ nj^Anj,B- 
{i,J) 



(2) 



where the sum over / runs over the superccUs on the lat- 
tice. Eq. ([2]) has the generic form of a two-band Hubbard 
model with the onsite intra-orbital Coulomb repulsion U 
and the inter-orbital Coulomb repulsion U' = V . The 
superlattice potential A^ plays the role of a crystal field 
and gives rise to a crystal field splitting A = A^ — A^. 
The last term with V' ~ V accounts for the inter-orbital 
Coulomb interaction between the neighboring supcrcells. 
Note that the U and U' are free parameters in our the- 
ory and not related by the Hund's rule coupling Jh in 
contrast to the usual multiorbital Hubbard modeli^. 



III. SUPERSTRUCTURE MOTT TRANSITION 
IN THE IONIC HUBBARD MODEL A / 0, V ^0 

We first focus on the ionic Hubbard model in the ab- 
sence of the inter-site Coulomb correlation {V = 0). The 
correlation effects will be treated by the Gutzwiller pro- 
jected wavefunctions. We first review briefly this vari- 
ational wavefunction approach to the Hubbard model 
and the Gutzwiller approximation as a semi-analytical 
method to carry out the projection. The virtual clus- 
ter Gutzwiller approximation will then be introduced to 
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FIG. 2: Dividing the original lattice into clusters. The dark 
dots represent the original lattice sites while the gray areas 
are the clusters, (a) a two-site supercell and (b) a four-site 
supercell. Here I and J label the position of each cluster. The 
dotted circles represent equivalent selection of the superceUs. 



treat the superlattice modulations induced by the ionic 
potential A^ ^ 0. 



Gutzviriller approximation to the Hubbard 
Model 



For the Hubbard model with only onsite repulsion U, 
the Hilbert space is a direct product of the single-site 
Fock states: one empty, two singly occupied (either by 
an up or a down spin), and one doubly occupied. The 
Gutzwiller projected wavefunction approach is a way 
of building correlation effects from the non-interacting 
wavefunctions by reducing the statistical weight of dou- 
ble occupatior>i^. The trial wavefunction can be written 
down formally as |^| ~ -P|*o), with |\l/o) being the non- 
interacting Slater determinant wavefunction and P the 
Gutzwiller projection operator. The projection opera- 
tor P ~ g^, where D = X^i^-iT^ii i^ ^^c double occu- 
pation operator of the system and 5 < 1 a variational 
parameter to be determined by minimizing the ground 
state energy. When minimizing the ground state energy 
{'i\H\'^)/{'^\'^), the average interaction energy is simply 
given by Ud, with d being the average density of double 
occupancy. Since there is a one to one correspondence 
between d and g, one can make d the variational param- 
eter instead of g. The task is then to derive an expres- 
sion for the kinetic energy in terms of d. Unfortunately, 
an analytical expression for this expectation value is not 
available other than in one and infinite dimensions'^. In 
two and three dimensions one has to resort to approx- 
imate methods. The Gutzwiller approximation relates 
the expectation values in the projected state with that of 
the uncorrelated state by a simple multiplicative factor, 
the Gutzwiller factor (GF). This is achieved by ignor- 
ing inter-site correlations in evaluating the expectation 
values with I 'J), 






*)/(*!*) 






c,.|*o), (3) 



where Qa is the GF for the kinetic hopping term and 
translational invariance of the lattice is assumed. The 
exact form of §„ depends on the trial wavefunction \^a) 
and can be obtained by considering how double occu- 
pation affects the hopping process. If \^q) is taken as 
the Slater determinant state without further symmetry 
breaking, g„ takes a simple form in terms of the uni- 
form particle density and the density of double occupa- 



tion^ii^, 



9a 



{ricj — d){\ + d — n) 



1/2 



d{ns - d) 



.1/2 



71^(1 - ric) 



(4) 



As a result, the minimization of the ground state en- 
ergy for the original Hubbard Hamiltonian translates to 
the minimization of the following rcnormalizcd mean field 
Hamiltonian. 



HGA^~Y.9lt 



V^ia^jo- ' 



'^X„{n,„-n^)+NUd, (5) 



where the hopping integrals are renormalized by a factor 
of g^ from the original ones. Here N is the number of 
sites and Act is a Lagrange multiplier that keeps the par- 
ticle density unchanged before and after the projection. 
Under the Gutzwiller approximation, the ground state 
energy of the original Hubbard model is approximated 
by that of the mean field Hamiltonian Hqa- The cor- 
relation effects are treated nonperturbatively since there 
are no self-energy corrections that scale with U in con- 
trast to the weak coupling Hartree-Fock approach. Since 
the bandwidth (hopping) is renormalized in Hqa, this 
approach is also known as the renormalized mean field 
approximation^^. 



B. The Virtual Cluster Gutzwiller Approximation 

Several authors have extended the original Gutzwiller 
approximation to handle lattice translational symme- 
try breaking using different schemes^ii^iS^i^. In 
the spatially unrestricted Gutzwiller approximation 
(SUGA)^i^, the GF for the hopping renormalization is 
generalized to depend on the local charge densities and 
double occupations at the sites connected by the hopping 
process. Specifically, gija for the hopping from site i to 
site j is a product gi^ ■ gj„, with gi„ given by 



[riic, - di) (1 + di - rii) 



1/2 



^i \J^i<7 



-,1/2 



(1 



ni„ (1 - n^o-) 

(6) 
where di is the double occupation at site i, riic is the 
spin dependent density, and n,; = ^^ riia. The latter has 
spatial modulations due to the presence of supcrlattice 
potential A^. The Mott transition in a given superstruc- 
ture at appropriate filling fractions can be analyzed using 
the SUGA and the GF in Eq. ®. 

In order to gain more physical insights into the prob- 
lem and establish the notion of superstructure Mottness, 
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FIG. 3: (Color online) The local Fock space of a two-site 
cluster. The expectation probability of each state and the 
corresponding interaction energy are also given. Here Et = 
U + 2V and Ef = 2U + 4V. 



we develop a virtual cluster approach. For a periodic 
supcrlattice potential A; ^ 0, the original lattice is di- 
vided into supercells. It is important to note that there 
are many equivalent ways to choose the supercell clus- 
ters that cover the whole lattice. A particular two-site 
supercell and a four-site supercell coverage of the original 
lattice is shown in Figs. [UJa) and[Ub). Thus, a lattice 
site is labeled by (/, i) where the capital letter / spec- 
ifies the supercell and i denotes the position inside the 
supercell. The Hamiltonian in Eq. ([1]) can be written as 
H = Hk + iJiocai, where the "kinetic" part Hk and the 
"local" interaction part i?iocai are given by, 

Hk=- Yl ^vC^HaCjja, (7a) 

local =U ^ nii^fiiii + ^ A^fiii. (7b) 

I.i lA 



H, 



In the spirit of Gutzwiller projected wavefunction, the 
effect of i/iocai is to reduce multi-occupation in each su- 
percell. For example, the Fock states of a two-site super- 
cell are shown in Fig. [Sj where 11 of the 16 states in a 
cluster have more than one particle. The operators for 
these multiplets are denoted by M^, 7 = 1, . . . , F, with F 
the dimension of the multi-occupancy space (F = 11 in 
the two-site cluster example). A trial wavefunction can 
then be written as j^f) = P|^o), where the projection 
operator P is a product of the cluster projection opera- 
tors, Pi, i.e., P = Y\j Pj. The cluster projection operator 
Pj reduces multi-occupation within each cluster / and is 
written down as. 



Pr = 



I[yt U 



(8) 



7=1 



where the supercell index / is omitted on the right hand 
side because of translation symmetry on the supcrlattice. 
Here the y's are fugacitics that keep charge densities un- 
changed upon projection at each site and the g-y^s are 



variationally determined weighting factors for the pro- 
jection. 



The hopping processes in Eq. (j7ap is then renormahzed 
by the projected wave function as in Eq. Q and the GF 
is calculated by ignoring inter-cell correlations. However, 
it is important to note that for a fixed choice of a super- 
cell coverage of the lattice, i^K includes hopping between 
different supercells/clusters (/ ^ J) as well as hopping 
within the same cluster (/ = J). Conventional cluster 
approaches, such as the one used by Lechermann et al., 
do fix the choice of the supercells. This, however, leads 
unavoidably to dimcrization since the inter-cluster hop- 
ping and intra-cluster hopping cannot be treated on equal 
footing and acquire different renormalization factors^S. 
Here we propose a virtual cluster approach. We do not 
fix a real space cluster coverage of the lattice in a pri- 
ori. Instead, the choice of the clusters is associated with 
the Hilbert space and arises when evaluating the expec- 
tation values of the hopping terms using the projected 
wavefunction in Eq. ([3]). Indeed, in order to implement 
the Gutzwiller approximation, the hopping term must 
be between different clusters (i.e. inter-cluster) and the 
correlation between the clusters must be ignored. Our 
insight is that for any two sites connected by the hop- 
ping process on the original lattice, one can always find 
a particular choice of superccU coverage such that the 
two sites reside in two different supercells. As an ex- 
ample, choosing the cluster marked by the dotted circles 
in Figs, dja) and ^h) turns the sites within the cluster 
(gray) into ones residing in different supercells. 



Specifically, consider the hopping terms in Eq. (j7ap . 
In the term where / = J, it is always possible to choose 
an equivalent cluster coverage such that the hopping is 
between sites in different clusters I' ^ J' . When evaluat- 
ing the expectation value of this hopping term using the 
projected wave function, one can simply rearrange the 
product of the cluster Fock states and write P = Y\j, Pi> 
and arrive at the same result as the case when I ^ J . 
It is in this sense, we term this approach as the vir- 
tual cluster Gutzwiller approximation (VCGA). One of 
the most notable advantages of this approach is that the 
same Gutzwiller renormalization factor is obtained for 
the hopping integral between any two sites regardless of 
the choice of the supercell coverage of the original lat- 
tice and, as a result, the difficulty associated with the 
unphysical dimcrization tendency is removed. 



The remaining steps for deriving the GF is the same 
as that of the multi-orbital Hubbard model with only 
on-site density-density interactions^^. In this sense, one 
can regard each supercell as an artificial "atom" and the 
lattice sites within the cluster as the associated atomic 
"orbitals" , and we will use orbitals and sites within a 



cluster interchangeably. We obtain 



jCjpa) = 



(9) 



{P!P])0 

(W,„,W7^,)o , t , 

(W/„^)0 {Wjpa)o i 

{Pfh {P])o '^' 

9laagjl3a{Cja„CJi3a)Q, 



where (■■•)o denotes expectation values in the unpro- 
jected state jvlfo) and the equahties hold only when 
inter-cluster correlations are ignored in the spirit of the 
Gutzwiller approximation. Here W/qct is defined as. 



rCjf3cr)Q 



J 
"laij 



Pl^ca 



Pi, 



(10) 



and giacr ~ {'^iaa)o/{Piaa)o- Bccausc of translational 
symmetry of the superlattice, giaa is independent of the 
supercell index /, i.e., gjaa- = gaa- 

In the rest of this section we focus on a two-site clus- 
ter on a square lattice, compatible with a superlattice 
potential 



-(-1) 



i:c+iv 



(11) 



The two sites inside a cluster, i.e. the two orbitals are la- 
beled as A and B. The A in the above equation has the 
meaning of a crystal field splitting between the orbitals. 
The multi-occupation projection operators M.^'s are as 

follows (suppressing the supercell index), Daa ~ DaEa 
- intraorbital doublon projection operator for the state 
with a doubly occupied a orbital (created by Da) and an 
empty orbital (created by Ea) ; -Dctct' = QaaQaa' - in- 
terorbital doublon projection operator for the state with 
two singly occupied orbitals of spin projections a (cre- 
ated by Qacr) and tr' (created by Qsa')] Taa = QaaDs 
~ triplon projection operator for the state with a singly 
occupied a orbital with spin a and a doubly occupied 
orbital; and F = D^Da - quadruplon projection oper- 
ator for the state with both orbitals doubly occupied. 
The projection operators for singly occupied states are 
denoted by Pacr where the a orbital with spin a is oc- 
cupied while the other orbital is empty, and project- 
ing to the configuration where both orbitals are empty 
is accomplished by £^ = EaEa- The expectation val- 
ues of these operators in the projected state will be de- 
noted by a lower case letter. For example, {Taa) = taa, 
see Fig. [3] The definition of the single-orbital projec- 
tion operators are the doubly occupied £)„ = fia^hai, 
singly occupied Qaa = ^acr(l — tt-uo) and empty orbital 
Ea = (1 — 77,o,|)(l — fia^). The expectation values of these 
single-orbital projection operators in the projected states 
are also denoted by a lower case letter, e.g., (Qaa) = <laa- 
It is seen from Eq. © that the GF, gaa, depends on the 
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FIG. 4: (Color online) Comparison of ground state energy of 
the Hubbard model at quarter filling treated with Gutzwiller 
approximation and variational Monte Carlo^. 1/ = A = 0, 
W = At. 



fugacities yaa and g-y. By expressing yaa and g-y in terms 
of the occupation probabilities, we obtain the Gutzwiller 
factors^, 



gAa = ^/epAa + \/dAAPAa + \J dcrcrPBa + \Jd^ 



PBo 



I \JnAcr{^ -riAcr), (12a) 



ffBcr = y^epBcr + y/dBBPBa + \/daaPAa + \/ dsaPAS 



+ \/ dAA^Ba + \J das^Aa + \J dsstAS + V^Bff/ 

l\/nBa{^-nBa)- (12b) 

In addition, the occupation probabilities satisfy the fol- 
lowing constraints, 

PAo = nAcj " dAA - 2J daai - t Aa - 2_j ^^'^' ^ •^' 



(13a) 



PBa = nsa 



— dsB — 2_^ '^o'''J ^ ^Ba — 2_^ ^Aa' — /, 



(13b) 



1 - ^ 71a + ^ daa + ^ d^a' + 2 ^ i^cr + 3/, 

(13c) 



where Eqs. (|13ap and (|13bp are just alternative expres- 
sions of the electron densities in terms of the occupation 
probabilities and Eq. p3cp is the completeness of the 
cluster Fock states. These results are consistent with 
those obtained for the multiorbital Hubbard modeli^. 



As an important check for the validity of the VCGA, 
we verified that the expressions in Eqs. (|12ap and (J12bp 
are equivalent to the results from the SUGA when inter- 
site correlations are ignored. To show this, one can make 
use of the fact that in the absence of inter-site corre- 
lations the multi-occupation probabilities in the virtual 
cluster approach arc simple products of the occupation 
probabilities in SUGA. For example, the above defined 
triple-occupation tAu is the product of the double occu- 
pation d-B and single occupation qAa-i i.e. tAa- = ^b ■lAa- 
By substituting all the multi-occupation probabilities in 
Eqs. (|12al and I12bp one recovers the results given in 
Eq. (El). 



C. Superstructure Mott Transition in the Ionic 
Hubbard Model 

We next present the results for the ionic Hubbard 
model obtained with the VCGA at quarter filling for the 
superlattice potential of Eq. (fTTjl . Since the system is far 
away from half filling, no magnetic instability is present 
and we will focus on the paramagnetic solutions. The 
effective Hamiltonian in the VCGA discussed above is 
given by, 

HgA = ^ gaag(3crS,al3{k)cl^^^Ck,pa + Y^ Aaflfc.a 
k.ajS k,a 

+ NcUC^daa +^tao- +2/) + ^ Xaa{nk,aa -riaa), 
a acr k,a<T 

(14) 

where Nc is the number of clusters and Acer's are the 
Lagrange multipliers that keep the particle density in 
each orbital unchanged following the projection. They 
originate from and have a one to one correspondence to 
the fugacities in the projection operator in Eq. ([8]). The 
ground state properties are determined by minimizing 
(Hga) and solving for ie,paa,daa,daa' ,taa, f) and A„<j 
sclf-consistently. 

To test the algorithm, we first compute the ground 
state energy per site at quarter-filling as a function of 
U/W, where M^ = 4Hs the half-bandwidth, for the Hub- 
bard model (i.e. set A = 0). Fig. |4] shows that the later 
agrees remarkably well with the results obtained using 
the variational Monte Carlo method to carry out the pro- 
jectioEk^. From the perspective of superstructure Mot- 
tness, it is clearly seen from Eq. (fT4)) that the U term 
alone docs not reduce the interorbital double occupancy 
daa' ■ As a consequence, U alone cannot drive a Mott 
transition away from half filling since the Gutzwiller fac- 
tors in Eqs. (|12ap and (J12bp are always finite and the 
system is in a correlated metallic phase. 

Switching on A ^ 0, we determined the phase dia- 
gram of the ionic Hubbard model shown in Fig. EJa) on 
the C/-A plane where a CDW metal at small {U, A) and a 
charge ordered insulator for large ([/, A) are separated by 
a continuous superstructure Mott metal-insulator tran- 
sition. The notion of superstructure Mottncss offers a 



physical explanation of the phase structure through the 
analogy to multi-orbital Mott-Hubbard systems. A finite 
crystal field splitting A induces orbital order ua ^ riB, 
which corresponds to charge ordering in the ionic Hub- 
bard model. However, below a critical strength of A, 
as seen in Fig. Efa), the system remains in a metallic 
charge density wave (CDW) state no matter how strong 
U is. This can be understood since the two bands derived 
from the two orbitals are not separated in this case and 
the inter-orbital hopping is not suppressed due to the ab- 
sence of an inter-orbital repulsion [/', which is related to 
the inter-site correlation V as shown in Eq. ^. 



A sufficiently large crystal field splitting A drives the 
system from a CDW metal to a charge ordered insulating 
state for large enough U. This happens because the two 
bands gradually separate with increasing orbital polariza- 
tion 5n = nB — nA- The lower band becomes half- filled 
since the system is at quarter-filling. A large enough 
U drives a Mott transition in the half-filled lower quasi- 
particle (QP) band. Indeed, as ua -^ and pAa -^ 0, 
Eq. (|13ap implies that all the double occupancies ap- 
proach zero. This leads to e — *■ at quarter filling from 
Eq. p3c[) . On the other hand, ub ^ such that from 
Eq. (|12bp the GutziwUer factor gBa -^ resulting in an 
orbital selective Mott transition shown in Fig. O^d) . In 
Figs.[5jb),[3c), the doublon densities and the inverse ef- 
fective mass of the lower QP band m/m* , which equals 
the QP spectral weight Z, are plotted as a function of A 
at a fixed large U = 3W. The continuous transition at 
Ac = 2AW is a Brinkman-Rice type transitional. The 
vanishing of the coherent QP spectral weight and the di- 
vergence of the effective mass (case shown in Fig. [Ub)) 
on approaching the transition from the metallic side are 
clear signatures of Mottness. 



The curvature of the phase boundary at large ionic 
potential A can also be understood. In this regime, the 
crystal field splitting is much larger than the hopping t 
and the bands are well separated. A simple calculation 
shows that the bandwidth w of the lower QP band de- 
creases with increasing A, w ^ t'^/A. As a result, the 
critical Uc for the Mott transition, which depends on the 
ratio of U/w^ decreases as can be seen in the phase dia- 
gram Fig. IHa). The evolution of the QP band dispersion 
as a function of A is shown in Fig. [DJd) at fixed large 
U = 3W. It clearly demonstrates that the Mott transi- 
tion in the ionic Hubbard model belongs to class (b) of 
the superstructure Mottness depicted in Fig. [TJ 



While the NN hopping t translates into interband hop- 
ping in the two-band model, the NNN hopping t' de- 
scribes intraband hopping. We find that including a 
t' = —0.2t moves the phase boundary toward smaller 
(A, [/), as can be seen in Fig. O^d) where a continuous 
transition takes place at a smaller Ac = 1.9VF. 
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FIG. 5: (Color online) (a) Phase diagram on tlie U-A plane 
for V = and t' = 0. Doublon densities (b) and inverse 
effective mass (c) as function of A for fixed U = 3W . (d) 
Band dispersion for different A, t' = — 0.2i. k = (7r,0) -^ 
(V2,^/2)^ (0,0)^(^,0). 



IV. SUPERSTRUCTURE MOTTNESS IN THE 
t-U-V MODEL 

In this section, we study the superstructure Mott tran- 
sition due to the inter-site Coulomb interaction V in the 
extended Hubbard, or the t-U-V model. To this end, the 
superlattice potential A is set to zero. Superstructure 
Mottness in the presence of both V and A will be stud- 
ied in the next section. The basic difficulty introduced 
by the inter-site interaction V is that the Hilbert space 
becomes nonlocal, making it inaccessible to the conven- 
tional Gutzwiller projected wavefunction. We show that 
the VCGA offers a natural nonperturbativc treatment 
of the Wigner-Mott transition within the framework of 
superstructure Mottness. 



Virtual cluster Gutzwiller approximation for 
the l/-term 



In the Gutzwiller approach, one would like to treat 
the interactions as a projection that reduces multi- 
occupation from an unprojected Slater-determinate 
state, i.e. j^f) = Pj^I'o). However, the projection oper- 
ator P cannot be written down in closed form, since the 
Hilbert space is infinitely connected by the inter-site V in 
Eq. 11]). Our strategy is to truncate the nonlocal Hilbert 
space as a product of that of the supercell clusters such 
that the U and V within a cluster serve as the intra- 
orbital U and inter-orbital U' just as in a multi-orbital 
Hubbard model, which can be treated nonperturbatively. 
The projection operator P can thus be written approxi- 
mately as the product of the cluster projection operators. 
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FIG. 6; (Color online) Comparison of the ground state energy 
obtained by VCGA and Hartree approximation of l/-term at 
U = W. 



Pi, i.e., P ~ Y\j Pj. This algorithm can be accomphshed 
by the VCGA for the ionic Hubbard model discussed in 
section III. A, with a few extensions. 

Although we do not fix a specific supcrcell structure in 
the virtual cluster approach, it is useful to write down the 
Hamiltonian for a given coverage of the lattice. Due to 
the non-local nature of the I^-term, the interaction part 
of the Hamiltonian can no longer be written as a direct 
summation over independent supercells. Let us write the 
Hamiltonian as the sum of two parts, similar to Eqs. ((7a|) 
and (l7bl). 



I A I ji'^J-j 



(15a) 
(15b) 



The second term in Eq. (|15bp contains both the intra- 
cluster as well as the inter-clustcr F-interactions. Con- 
sider now the projection of the hopping term. In the 
VCGA discussed before, for a given hopping term be- 
tween two sites, one can always choose a cluster coverage 
such that the two sites reside in two different clusters. 
The expectation value of the hopping term between the 
projected states can be evaluated according to Eq. ©. 
Moreover, to be consistent with the Gutzwiller approx- 
imation, the inter-cluster correlation should be ignored, 
which means that the inter-clustcr V in Eq. (|15bp should 
be turned off accordingly. The form of the cluster pro- 
jection operator is therefore formally identical to that 
in Eq. ([8]) for the ionic Hubbard model, since the intra- 
cluster V-tcim in Eq. (|15bp only leads to different ener- 
gies of the multiplets in the cluster Hilbert space. An 
explicit example of the F-dependent cluster Fock state 
energies is given in Fig. [3] for the case of a two-site clus- 
ter. Consequently, the Gutzwiller hopping renormaliza- 
tion factors for all hopping amplitudes have the same 



form and arc given by Eqs. (|12ap and (J12bp in terms 
of occupation probabilities for two-site clusters. The 
inter-clustcr I^-term will take the same value since the 
bond becomes intra-cell under a different choice of the 
virtual cluster configuration. In short, the advantage 
of the virtual cluster formulation is that all bonds con- 
nected by hopping and the inter-site Coulomb interaction 
V are treated on equal footing and remain equivalent af- 
ter Gutzwiller projection without artificially breaking the 
lattice translation symmetry. 



B. Wigner-Mott transition as superstructure 
Mottness 



We now present our results obtained for the extended 
Hubbard t-U-V model at quarter-filling using the VCGA 
with two-site clusters. Following the above discussion, 
the renormalized mean-field Hamiltonian in the VCGA 
is given by, 

k,af3 

+ N.UiY, ^"" + !]*«- + 2/) 



aa' OiO 

k,a(T 



(16) 



where the Gutzwiller factors are given by Eqs. (|12ap and 
(J12b[) with the occupation probabilities satisfying the 
fermion counting and completeness equations (|13aj|13c() . 
The paramagnetic ground state is obtained by minimize 
(Hqa) with respect to the variational parameters just as 
in the ionic Hubbard model studied in Section III. 



1. Comparison to Hartree approximation 

We first demonstrate that our strong coupling, non- 
perturbativc VCGA treatment of the inter-site correla- 
tion V is fundamentally superior to the weak-coupling 
Hartree approximation of the y-term with a single-site 
Gutzwiller projection of the on-site repulsion U. To this 
end, we compare the ground state energy of the quarter- 
filled t-U-V model at U = W. Fig. [6] shows that while 
the ground state energy obtained by VCGA is consistent 
with the weak coupling Hartree result for small values 
of V, which is in fact a nontrivial check for the validity 
of the VCGA, it is significantly lower for all values of V 
larger than about 10% of the bandwidth. The ground 
state energy approaches zero for V ^ W in the Hartree 
approach, signifying a transition to the charge ordered 
insulating state. The result of the VCGA clearly shows 
that the Hartree approximation grossly overestimated 
the charge ordering tendency. This is because the weak- 
coupling Hartrec-decoupling gives rise to a self-energy for 
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FIG. 7: (Color online) (a) Phase diagram of the t-U-V model 
on the U-V plane, (b) Phase diagram in the F-Hartree ap- 
proximation. Doublon densities (c) and inverse effective mass 
(d) as functions of V at fixed U — 3W. 
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FIG. 8: (Color online) (a) Phase Diagram of the t-U-V-A 
model on the V-A plane at fixed U = 3W. The circle on the 
boundary separates discontinuous and continuous transitions, 
(b) Charge difference Sn = ub — ua as a function of V with 
A = i and U = 3W^. (c) The variational energy as a function 
of 5n at A = i, 1/ = 2.6i, and U = 3W. 



the fcrmions that scales with V and becomes unphysical 
when V is large. On the other hand, in the nonperturba- 
tive VCGA, there is no self-energy cost that scales with 
V . Thus the uniform metallic ground state turns out to 
be much more robust against inter-site correlations. 

While it is clear that the results of the VCGA can 
be systematically improved with increasing cluster size 
to treat the inter-site correlations even better, we show 
below that it already captures the essence of the super- 
structure Mott metal-insulator transition with two-site 
supercells. 



2. Wigner-Mott transition in the t-U-V model 

The phase diagram obtained using the VCGA at 
quarter-filling is shown in Fig. [7{a) on the U-V plane. 
It contains a uniform strongly correlated metal at small 
([/, y) and a Mott insulator at large {U^V). Since the 
crystal field splitting is absent, the two degenerate QP 
bands narrow and undergo a continuous Mott transi- 
tion simultaneously upon approaching the phase bound- 
ary from the metallic side, corresponding to the first 
case of the superstructure Mottness transition shown in 
Fig. [TJa). Figs.[7l^c) and[7jd) show that double occupa- 
tion, the (inverse) QP mass renormalization m/m* and 
the QP coherence weight decrease continuously with in- 
creasing V , which are clear signatures of Mottness. and 
vanish at the critical K = 1-9VF for U = iW . The 
strongly correlated metallic phase is stable against the 
two-sublattice antiferromagnetic (AF) spin density wave 
order because the quarter-filled Fermi surface is away 
from the AF zone boundary. It is remarkable that for 
degenerate bands, the Mott transition driven by V takes 
place with uniform charge density without forming a 
Wigner lattice, suggesting that superstructure Mottness 



is a more suitable description than a Wigner crystal. 
Since the Mott insulating state has the same energy as 
the (0, 1) charge ordered classical Wigner crystal, an ar- 
bitrarily small crystal field (i.e. superlattice potential) 
would trigger a first order transition to an orbital or- 
dered state, corresponding to a charge ordered state in 
the original model. Interestingly, if the two bands have 
different bandwidths, corresponding to anisotropic NNN 
hopping t'j^ 7^ ^Bi the Mott insulating state emerges at 
large U and V continuously from a CDW metal with a 
finite density polarization 5n. 

It is useful to further compare these results to previ- 
ous studies where the l^-term is decoupled by Hartree 
approximationi^ii^ii^. The phase diagram in the latter 
case obtained by the single-site Gutzwiller approxima- 
tion is shown in Fig. [7Jb), where a uniform metal at 
small V transforms into a CDW metal and eventually 
to a charge ordered insulator at large V . The topology 
of this phase diagram is consistent with that obtained 
by the dynamical mean field theory at finite tempera- 
tures with a semi-circular density of states^^. However, 
the CDW metal phase is absent in our VCGA approach 
that treats U and V on equal footing. It is likely that in 
the Hartree approximation, the fermion self-energy scal- 
ing with V overemphasizes the symmetry breaking CDW 
ordering tendency at large V as discussed above. 



V. SUPERSTRUCTURE MOTTNESS WITH 
COMBINED A AND V 

Finally, we turn to the most general case of the t-U- 
V-/S. model and study the interplay between the crystal 
field splitting A and the extended Coulomb interaction 
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V. Treating the V-tcTin using the VCGA in the same 
way as in Sec. IIVI and putting the A-term in the lo- 
cal part of the Hamiltonian in Eq. (J15bp . we arrive at 
the identical form of the Gutzwiller renormalization fac- 
tors. The only difference is that the crystal field term 
should be added in the renormalized mean-field Hamil- 
tonian in Eq. P^ . Fig. [SJ^a) displays the phase diagram 
at U = 3W for the quarter-filled t-U-V-A model. In 
contrast to the V = case shown in Fig. [Slja), the co- 
herent motion in the strongly correlated metallic state 
is further suppressed by the inter-site correlation V such 
that a small (infinitesimal) A can drive the system to the 
charge ordered insulating state, provided V is large. Sim- 
ilarly, the superlattice potential A is seen to enhance the 
correlation effects such that the localization and local mo- 
ment formation can be induced by a moderate strength 
of the inter-site Coulomb interaction V . 

The phase boundary between the CDW metal and the 
charge ordered insulator in Fig. [5Ja) changes from a dis- 
continuous Wigner-Mott transition at small A to a con- 
tinuous Brinkman-Rice transition above a critical A. The 
existence of two types of metal-insulator transition when 
both V and A are present is consistent with a recent 
DMFT study of a two-band Hubbard model with U = U' 
and A 7^ 0^. Fig. [H{b) shows the charge density differ- 
ence as a function of V at A = t, displaying a discon- 
tinuous jump at a critical Vc = 2.6f. The discontinuous 
transition is a consequence of degenerate minima in the 
ground state energy for fixed {V, A). In Fig.[5Jc), we plot 
the variational ground state energy -Evar as a function of 
the density difference Sn at V ~ 2.6t and A = i close to 
the phase boundary. The minimum at dn = 0.67 corre- 
sponding to a CDW metal is nearly degenerate to that 
of the charge ordered insulator at Sn = 1. Increasing V 
and A further will make the energy of the insulating state 



the global minimum and trigger the first order transition. 
The existence of two types of localization transition into 
the local moment phase is an important and unique prop- 
erty associated with superstructure Mottness. 



VI. CONCLUSIONS 

We have shown that the superlattice potential in an in- 
homogeneous electronic state and the inter-site Coulomb 
repulsion increase the correlation effect and the tendency 
toward localization and local moment formation. Map- 
ping to the generic multiband Hubbard model revealed 
a deeper connection among the strongly correlated inho- 
mogcneous electronic states, the Wigncr-Mott physics, 
and the multiorbital Mott-Hubbard physics under the 
unified notion of superstructure Mottness. An unbiased 
treatment of both the on-site and inter-site Coulomb in- 
teractions, e.g. the virtual cluster approximation devel- 
oped here within the framework of the Gutzwiller ap- 
proximation, is essential as demonstrated for the quarter 
filled Hubbard model. We expect that doping away from 
such unconventional insulating states will further reveal 
behaviors of superstructure Mottness, such as the coex- 
istence of local moment and itinerant carriers. At dilute 
electron densities, the competition between itinerant fer- 
romagnetism and local moment formation may also arise 
due to superstructure Mottness in inhomogeneous elec- 
tron systems. 
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